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We investigate the possibility to detect the integrated Sachs-Wolfe (ISW) effect by cross- 
correlating 21-cm surveys at high redshifts with galaxies, in a way similar to the usual CMB-galaxy 
cross-correlation. The high-redshift 21-cm signal is dominated by CMB photons that travel freely 
without interacting with the intervening matter, and hence its late-time ISW signature should cor¬ 
relate extremely well with that of the CMB at its peak frequencies. Using the 21-cm temperature 
brightness instead of the CMB would thus be a further check of the detection of the ISW effect, 
measured by different instruments at different frequencies and suffering from different systematics. 
We also study the ISW effect on the photons that are scattered by HI clouds. We show that a 
detection of the unscattered photons is achievable with planned radio arrays, while one using scat¬ 
tered photons will require advanced radio interferometers, either an extended version of the planned 
Square Kilometre Array or futuristic experiments such as a lunar radio array. 


The integrated Sachs-Wolfe (ISW) effect pj{3] is a 
gravitational redshift due to the evolution of the gravita¬ 
tional potential as photons pass through matter under¬ 
and over-densities in their path from the last scatter¬ 
ing surface to us. In an Einstein-de Sitter universe, the 
blueshift of a photon falling into a well is cancelled by 
the redshift as it climbs out, giving a zero effect. In 
the presence of a time variation of the local gravitational 
potential, due to a dark energy component or a devia¬ 
tion from Einstein’s theory of gravity, however, potential 
wells are modified and photons will experience a blue- or 
red-shift, leading to a net change in photon temperature 
which accumulates along the photon path. This effect 
translates into temperature anisotropies proportional to 
the variation of the gravitational potentials. The ISW 
effect has been detected mn] through cross-correlation 
of CMB maps at GHz-frequencies with galaxy surveys 
and used to constrain cosmological parameters [IMS]. 

In a similar way, the “21-cm background” at high- 
redshift, which is sourced by CMB photons in the 
Rayleigh-Jeans tail of the blackbody spectrum that travel 
through neutral hydrogen clouds, will experience an ISW 
effect from the evolution of gravitational potential wells 
(see Fig. [I]). Stripped of its unscattered CMB contri¬ 
bution, the 21-cm signal represents a unique source of 
cosmological information from the epoch of reionization 
(EOR) and the dark ages 16. jTT] ; measuring at differ¬ 
ent frequencies, the large redshift volume of these two 
epochs can be used to probe the spatial distribution of 
neutral hydrogen p~8U20l . The 3D power spectrum and 
imaging of the 21-cm brightness temperature fluctuations 
can also be used for various other cosmological mea¬ 
surements (see e.g. [21] [22]). Given that mapping the 
21-cm signal is the target of several ongoing and near¬ 
future ground-based experiments such as the SKA 23] , 
as well as more futuristic experiments on the far side 
of the Moon [24, 25 or space-based radio arrays like 


RadioAstron [26], DARE [27], DALI [28], OLFAR [29] . 
DARIS [30], SURO [3T] , PARIS [32], that promise to 
reach very high accuracy, it is worthwhile to investigate 
how to maximize the scientific output of such measure¬ 
ments. In this Letter we briefly present the formalism 
for studying the ISW effect on the 21-cm signal and 
the experimental challenges for its detection, using the 



FIG. 1. Illustration: The radiative transfer of CMB photons 
through neutral hydrogen gas clouds induces fluctuations at 
21-cm frequencies (due to absorption or emission, depending 
on the relative temperatures of the inter-galactic medium and 
the CMB) which are the target of many ongoing and planned 
experiments. The majority of the signal, however, is com¬ 
prised of unscattered CMB photons at the Rayleigh-Jeans tail 
of its blackbody spectrum. All these photons later undergo 
line-of-sight blue- or red- shifting as they travel through evolv¬ 
ing gravitational potential wells (e.g. during A-domination). 
This integrated Sachs-Wolfe contribution to the fluctuations 
at 21-cm frequencies can then be efficiently extracted by cross- 
correlating the maps with low-redshift galaxy surveys. 
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cross-correlation of low-redshift tracers of matter (galax¬ 
ies) with maps of brightness temperature fluctuations at 
21-cm frequencies corresponding to high redshifts (see il¬ 
lustration in Fig. [l]). 

The dominant signal in these maps is of unscattered 
CMB photons, and therefore its late-time ISW signature 
will be very highly correlated with the signature at the 
peak CMB frequencies mentioned above. A complement¬ 
ing measurement at 21-cm frequencies would be promis¬ 
ing as it represents an independent detection of the ISW 
effect, measured with different instruments and contami¬ 
nated by different foregrounds. As the 21-cm background 
is set to be observed across a vast redshift range by up¬ 
coming experiments, we show that there should be ample 
signal-to-noise for this detection. There is, however, an¬ 
other feature that could be sought after, which is the ISW 
effect on those CMB photons that do interact with the 
neutral hydrogen clouds at high redshifts. Assuming the 
CMB fluctuations are efficiently subtracted from the 21- 
cm maps, we demonstrate that this signal can potentially 
be detected in the data as well, although most-likely re¬ 
quiring futuristic experiments. 

We parameterize potential perturbations in the confor¬ 
mal Newtonian (cN) gauge 

ds 2 = a 2 (r]) [—(1 + 2+ (1 — 2<F) Sijdx l dx j ] . (1) 

The (fractional) shift in the observed photon frequency 
is given by 

^ = (*. - *„) + d X (* + $). (2) 

The right hand side is frequency-independent. The first 
term, known as the Sachs-Wolfe term, is the potential dif¬ 
ference between the source location s and the observer’s 
location Q . The second term is an ISW line-of-sight in¬ 
tegral from the observer to the source, with y being the 
comoving radial distance. In this work we focus on the 
second term. There is also a second order effect, called 
Rees-Sciama [33], due to the evolution of clustering while 
the photon is traversing the potential well, that happens 
even in the absence of dark energy, but this is typically 
orders of magnitude smaller than the ISW effect, so we 
will neglect it. 

The observed 21-cm signature at frequency is will cor¬ 
respond to frequency v — Sis without the ISW effect. By 
conservation of photon number, the observed occupation 
number (per polarization state) f(y) would be 


/(0 = fo(y - Su) = f 0 (u) 


din f Sis 
dims is 


(3) 


where fo(is) is the occupation number without ISW. The 
fractional change is therefore given by 

<5/ d\nf 8v 

f d\nv v' U 


At a given observed frequency z/, the specific intensity 
I(y) (for both polarization states) is related to the oc¬ 
cupation number by I{is)dis = 2f(is)is 2 dis. We therefore 
find a fractional change (note that v is fixed) 


8I_ 

T 


sf 

f 


din f Sis 
dims is 


(5) 


The 21-cm signature is sourced by photons from the 
CMB background that traverse neutral hydrogen regions 
before reaching the observer. If interaction with neu¬ 
tral hydrogen clouds is neglected, the occupation number 
is given by a blackbody spectrum, but in the Rayleigh- 
Jeans regime 
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characterized by a thermodynamic temperature Tr 0 = 
2.725 K (at present redshift). It then follows that 
—din f/dims = 1, and therefore 
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Since we are in the Rayleigh-Jeans limit, the specific 
intensity is proportional to the temperature, and we 
reproduce the ISW contribution to CMB temperature 
anisotropies 
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If the interaction between the CMB and neutral hydro¬ 
gen is accounted for, the occupation number f(is ) will be 
distorted 
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where the mean brightness temperature T 2 i is given by 
T s (z)-T r (z) 
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( 10 ) 


corresponding to the redshift z = z(is) of CMB-atom 
interaction. Here r 2 i is the optical depth for the hyperfine 
transition [34] 
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Thus, the intensity fluctuation is given by 
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Eq. (12) is the main equation of this paper, and includes 
the two signatures we are targeting. The second term in 
the parentheses is suppressed by T 2 i/Tro nr 2 -nr 1 , 
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and is therefore generically much smaller than unity. As 
a first approximation, we can choose to neglect it and 
consider first the detection of the ISW effect on the un¬ 
scattered CMB photons, represented by the first term. 
Then, assuming the unscattered CMB fluctuations can 
be removed using external CMB datasets, we can es¬ 
timate if the ISW effect on the scattered photons can 
be detected as well. For this purpose, the ideal redshift 
would be z < 30, where the derivative SlnXA/dlnr' can 
turn out to be large (see e.g. Figs. 4, 7 of |T7]). This 
redshift corresponds to the very end of the dark ages 
and might be within reach of the SKA. The exact value 
of the second-term coefficient is highly model-dependent, 
though (mainly dictated by the X-ray and Lyman-<a emis- 
sivity), and we choose a fiducial (optimistic) value of 0.3 
for the purpose of examining its detectability. Since pho¬ 
tons from z ~ 1100 pass through an epoch with some 
residual radiation domination, they undergo an addi¬ 
tional early ISW redshift not experienced by those from 
z ~ 30. There will therefore perhaps be subtle differences 
between the temperature fluctuation in the scattered and 
unscattered photons. We leave a more detailed discussion 
of this effect to future work. 

The ISW effect can contribute significantly to the tem¬ 
perature fluctuations on large angular scales, but it en¬ 
hances only the low multipoles, and it is generally smaller 
than the temperature power spectrum, both at the CMB 
peak frequencies and at those corresponding to high red¬ 
shift 21-cm. Typically the ISW signal is weak and cross¬ 
correlation with tracers of the density field is necessary 
to bring it to measurable levels. 

We can write the cross-correlation power spectrum be¬ 
tween the surface density fluctuations of galaxies and 
temperature fluctuations as (for more details, see e.g. [2j 

13ISI ESI) 

cf = <a? ro O = 4tt J yA \k)Wf{k)Wf{k ), (13) 

where Wf and Wj are the galaxy and temperature win¬ 
dow functions respectively, and A 2 (k) is the logarithmic 
matter power spectrum today. 

The galaxy window function can be written as 

Wf{k) = J dz ^ b(z)D(z)j e [kx(z )\, (14) 

where ( dN/dz)dz is the mean number of sources per 
steradian with redshift 2 ) within dz , brighter than the 
flux or magnitude limit; b(z) is the bias factor relating 
the source to the mass overdensity, D(z) is the linear 
growth factor of mass fluctuations; ji is the spherical 
Bessel function of order t, and x( z ) is the comoving dis¬ 
tance. We use a redshift distribution and bias for galaxies 
as predicted for the SKA survey, but our results are not 
heavily dependent on them. Selecting a survey with the 
appropriate median redshift and the bias of the objects 


observed, one can in principle optimize the detection of 
the cross-correlation; a detailed study of this is beyond 
the scope of this paper. 

The window function for the ISW effect is 

wj{k) = I dx T x (z)W( X (z), X „)je[k X (z)}, (15) 

where \v corresponds to the central frequency targeted 
by the experiment. 

The errors for the cross-correlation can be computed 
as (see e.g. [35] ) 


°c* T = 


N 


(Cf) + [(Cf fl + n 9 1 ) + e)] 

( 2 £ + l)/ sky 


(16) 


where f s k y is the sky coverage of the survey, n g is the av¬ 
erage number density of sources, and 6 is the instrument 
noise. In the case of CMB photons this corresponds to 
the beam noise and it’s usually considered to be negli¬ 
gible for modern experiments; for 21-cm experiments we 
have [2T] 


2 7t/3(z) 
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where Av is the frequency bandwidth, /3(z) = 

( 27 r) 2 T s 2 ys (z)(^ ov )-i, with T sys « lSO^/lSOMHz)- 2 - 6 K, 
and £ cov = 2 ttV/X , where V is the length of the inter¬ 
ferometer baseline. Finally, we take Cj T to be the regu¬ 
lar CMB power spectrum when considering the signal-to- 
noise for detection of the ISW effect on the unscattered 
CMB photons and replace it with the 21cm power spec¬ 
trum when considering the effect on the scattered ones 
(assuming perfect CMB subtraction). 

In order to detect the ISW effect we need to cross¬ 
correlate temperature maps with galaxy catalogs. In the 
case of the CMB-galaxy correlation, 5 Cj T and hence 
the signal-to-noise ratio (SNR) of the cross-correlation 
depends, as in the “standard” ISW case, mostly on the 
sky coverage of the survey and the number of sources in 
the galaxy catalog. 

The error budget in the case of the scattered photons is 
dominated by the e term of Eq. (16), so the detectability 
of this cross-correlation largely depends on being able to 
reduce it. We study the SNR of the 21-cm-galaxy corre¬ 
lation as a function of the interferometer specifications. 

In Figure |2] we show the SNR for the unscattered (up¬ 
per panel) and scattered (lower panel) cases, where we 
used the predicted catalog of sources to be detected with 
the SKA. 


In the first case we can see that forthcoming and 
planned surveys such as ASKAP and the SKA could 
achieve a SNR of 5-6 (for specifications of these sur¬ 
veys and predictions about detecting the ISW by cross- 
correlating their sources with CMB experiments, see m 
l36U38] h Using multiple redshift bins, this signal-to-noise 
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could be further increased by a factor of a few. In the 
upper panel of Figure [2] we show the SNR for the cor¬ 
relation of LSS with unscattered photons, as a function 
of the area of the sky surveyed and the galaxy number 
density. 

In the latter case, a detection is, as expected, more 
difficult to obtain. The 5 term depends on a series of 
parameters of the 21-cm detecting instrument, such as 
the observing time, the frequency bandwidth, the frac¬ 
tional area coverage and the length of the baseline. In 
the lower panel of Figure [2] we plot the SNR for the 
cross-correlation of LSS with 21cm, scattered, photons, 
as a function of the fractional area coverage and baseline 
length. To calculate the 21-cm power spectrum, Cf 121 , 
we used CAMB Sources [39] with the appropriate fre¬ 
quencies and bandwidths. The results weakly depend on 
the galaxy survey used. Different surveys give slightly 
different results, but we do not find a dramatic change in 
the overall signal-to-noise ratio; targeting specific redshift 
ranges and objects could help. Other ways to improve the 
results could come from the use of a tomographic analysis 
in the galaxy catalog as well, and from the combination 
of different surveys (see e.g. [6, 14]). The combination of 
different surveys, in particular, has already been shown to 
improve the detection of the signal in the case of the LSS- 
CMB correlation. These optimizations could improve the 
ability of detecting the ISW from the LSS-21cm correla¬ 
tion, but a detailed analysis of these prospects is beyond 
the scope of this work and will be studied elsewhere. 

To conclude, we have investigated the possibility to de¬ 
tect the ISW effect using high-redshift 21-cm maps mea¬ 
surements. Our results show that currently planned in¬ 
struments will be able to detect the ISW effect on the 
CMB photons that pass freely through the neutral hy¬ 
drogen with a fairly decent sensitivity, allowing a cross¬ 
check to those measurements which use the CMB at peak 
frequencies. As for the effect on the scattered photons, 
currently planned experiments will not be able to pro¬ 
vide a clear detection of the signal, even via the cross¬ 
correlation with tracers of the underlying matter distri¬ 
bution such as galaxies. Future experiments, such as an 
improved SKA or a lunar radio array on the far side of 
the moon, however, could possibly grant a signal-to-noise 
that could allow a clear detection of this cross-correlation. 
Both measurements could be used as further tests of the 
standard cosmological model. 
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FIG. 2. Upper Panel : Signal-to-noise ratio for the correlation 
of LSS with unscattered photons, as a function of the area 
of the sky surveyed and the galaxy number density. Lower 
Panel : Signal-to-noise ratio for the LSS-21cm correlation, as 
a function of the fractional area coverage and baseline length. 
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